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1. 

Thin  Unsteady  rleavY   Jets 
Joseph  3.  I-Ie?.ler  and  riortimer  'veltz 

1 .   Introduction 

The  hydrodynamic  theory  of  jets  attempts  to  determine 
the  surfaces  of  a  jet  and  the  detailed  velocity  distribution 
within  a  jet.   Tlie  main  aifficulty  stems  from  the  nonlinear 
boundary  conditions  to  be  satisfied  on  the  unknown  jet 
surfaces.   >/hen  the  flow  is  restricted  to  be  two-dimensional 
and  steady,  and  surface  tension  is  neglected,  this  difficulty 
can  be  overcome  by  the  method  of  conformal  mapping  and  the 
surfaces  and  flow  can  be  found  exactly.   This  determination 
is  based  on  the  exact  equations  of  hydrodynamics,  om.itting 
viscosity. 

In  order  to  treat  more  general  problems  (unsteady, 
three-dimensional,  or  including  gravity)  the  simpler  hydraulic 
theor7y"  is  usually  em-oloyed.   In  this  theory  one  gives  up  the 
attempt  to  find  the  detailed  velocity  distribution  within  the 
jet,  but  instead  assumes  that  the  velocity  is  constant  on 
each  cross-section  of  the  jet.   Since  this  assumption  is 
incompatible  with  the  exact  equations  of  hydrodynamics,  it 
is  necessary  to  use  different,  approximate,  equations  based 
on  the  conservation  lax-fs  of  mechanics.   These  equations  can 
be  solved  and  yield  the  shape,  thickness  and  speed  of  the  jet, 
approximately.   (This  method  can  also  be  extended  to  take 
account  of  surface  tension,  as  will  be  shown  belov;.) 

Two  questions  which  immediately  arise  arc:   "i-rnat  is  the 
relationship  between  the  two  theories?"  and  ''How  can  the 
results  of  the  hydraulic  theory  be  improved?"   In  this  paper 
we  ansx-rer  these  questions  'oj   presenting  a  m.ethod  of  solution 
of  the  hydrodynamic  problem  as  a  series  in  powers  of  the  jet 
thickness  divided  by  som.e  other  typical  length  of  the  jet. 
The  first  torn  in  this  solution  is  found  to  be  the  solution 
given  by  the  hydraulic  theory,  thus  answering  the  first 
question.   The  higher  order  terms  in  the  series  yield 
corrections  to  the  hydraulic  theory,  thus  answering  the  second 
question. 
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We  vrere  led  to  the  solution  by  notins;  tlif-it  in  the 
problem  of  water  waves,  which  is  mathematically  sirailar  to 
that  of  a  jet,  both  gravity  and  time  dependence  are  taken 
into  account.   Because  of  the  similarity  of  the  two  problems, 
we  expected  that  the  methods  used  in  treating  water  waves 
could  also  be  applied  to  jets.   VJg  found  that  this  was  indeed 
the  case.   In  fact  precisely  the  same  kind  of  problem  arises 
in  \-jater   wave  theory  in  showing  the  connection  between  the 
exact  theory  and  the  so-called  non-linear  shallow  water  theory. 
(The  shallow  water  theory  is  the  hydraulic  theory  of  water 
Xijavo  s . ) 

The  mxthod  vjhlch  wo  cjiploy  involves  the  introduction 
of  dimensionless  variables  in  such  a  way  that  horizontal  and 
vertical  distances  are  stretched  'o'j   differezit  amounts,  just 
as  in  the  treatment  of  waves  in  shallow  water,   rlowever,  we 
must  take  account  of  the  fact  that  the  flow  is  not  essentially 
horizontal,  as  in  the  shallow  water  case,  but  is  parallel  to 
some  ujiknown  curve.   (In  flow  over  a  spillway,  which  vjo  also 
consider,  this  curve  is  the  bottom  and  is  therefore  known.) 
This  difference  is  taken  account  of  in  the  choice  of 
dimensionless  variables.   Once  appropriate  variables  have 
been  introduced,  the  solution  is  expressed  as  a  series  in  the 
ratio  of  the  horizontal  to  the  vertical  scale  factoro*   The- 
successive  terms  in  the  series  are  then  obtained  in  a 
straightforward  manner. 

This  procedure  is  carried  out  for  two  dimensional  jets. 
The  results  apply  to  both  steady  and  unsteady  jets,  v;ith  or 
vjithout  gravity.   Thus  the  initial  width,  speed  and  direction 
of  the  jet  may  vary  in  time.   In  the  resulting  zero- order  or 
hydraulic  theory  the  individual  "particles''  of  fluid  move 
along  parabolic  trajectories  as  if  they  were  freely  falling, 
and  the  pressure  throughout  the  jet  is  the  saie  as  the 
external  pressure.   In  the  case  of  flow  along  a  spillvray  or 
hill,  the  flov;  may  detach  from  the  bottom  if  the  bottom 
curvature  is  too  rrreat. 
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In  section  2  the  hydrodynamic  problem  is  formulated, 
and  the  diriensionlc  ss  variables  and  series  expansion  are 
introduced.   In  section  3  the  zero  order  solution  for  both 
steady  and  unsteady  jets  is  given,  and  found  to  agree  with 
the  hydraulic  theory.   In  section  l\.   the  first  order  solution 
for  steady  jets  is  obtained,  which  represents  an  improvement 
over  the  hydraulic  theory.   In  section  5  the  steady  and 
unsteady  flow  over  a  rigid  surface,  such  as  a  hill  or  spillway, 
is  treated.   In  section  6  the  hydrodynamic  formulation  is 
modified  to  include  surface  tension,  and  the  resu.lting  zero 
order  equations  are  foujid  to  coincide  with  those  of  a 
hydraulic  theory  which  includes  surface  tension.   The  zero 
order  solution  for  a  steady  jet  is  also  given,   (Appendix  I 
treats  vertical  stead;/  jets,  which  require  a  separate  analysis, 
and  both  zero  and  first  order  solutions  are  obtained.   A  very 
slight  modification  of  this  method  enables  us  to  treat 
vertical  jets  of  circular  cross  section,  but  this  is  not  given 
here . ) 

2 •   formulation  of  the  Problem 

We  consider  the  two-dimensional,  irrotational  flov;  of 
an  incompressible,  inviscid  fluid  bounded  by  two  free  stream- 
lines y  =  y.(x,'t)j  i  =  1,2.   This  flovj  may  be  called  an 
unsteady  jot.   The  coordinate  system  is  so  chosen  that  gravity 
acts  in  the  negative  y  direction.   We  assume  that  the  pressure 
has  the  same  constant  value  on  both  free  stream.lincs,  and 
this  constant  will  be  taken  to  be  zero,  without  loss  of 
generality.   The  horizontal  velocity  u,  the  vertical  velocity 
v,  and  the  pressure  p  satisfy  the  following  differential 
cqu.ations  and  boundary  conditions  [1]  (n  is  the  density  and 
g  the  acceleration  of  gravity): 
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u   +  V   =0 

X   y 


u_  +  uu_  +  vu_  =  -  p_/p 


(2.1)       V   +  uv_  +  vv^  =  -  p_/p  -  g 
t     X     y      y 

v_  -  {!_  =  0 


p  =  0  at  y  =  f^{x,t)  (i  =  1,2) 


y  _  +  y  „  ^  =  ^  ^■'^  y  =  y.-(^,t)    (i  =  1,2) 

it    Ix 


In  addition  to  those  Gqviations,  the  values  of  u,  v,  d,  y. 
must  bo  flven  at  -?ofne  cross- section  of  the  jot  (the  "nozzle"). 
These  conditions  v;ill  bo  considered  after  we  obtain  a  general 
expression  for  the  solution. 

We  nov/  intro'luce  dimensionloss  variables,  choosing  two 
constants  a  and  b  with  dimensions  of  length,  as  scale  factors 
in  the  horizontal  and  vertical  directions  respectively. 
These  constants  may  be  typical  dimensions  of  the  jet  in  the 
respective  directions.   The  dimensionless  variables  arc 
defined  by: 

x   =  ax  t={—  tp  =  agpp 

(2.2)  y   =   ay^(x,t)    +   by  u   =   J~  u 

y^(x)    =   ay^(x,t)    +   by^(K,t)  v   =  Jag   v 

In   these   equations  y  =   ay   (x,t)    is   the   equation   of   the   as   yet 

undetermined    center   line    of  the    jet,    i.u.      y    (x,t)    is    to  be 

chosen   such  that  v,  +  y^   =  2aY    .      From  this    it   follows   that 

"1      ■'  z  "^  o 
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(2.3)  yi(-,t)  =  -  y2(x,t) 

Inserting  the  new  variables,  equations  (2.2),  into 
equations  (2.1)  wo  obtain  (lotting  —  =  CC   ) 

QTu   -  V  u   +  V   =0 
X    -  ox  y    y 

<Tu .     -   7    ,  u      +  cTuu      -   Y     uu      +   vu      =   -  (Tp      +   v      p 
^    t        -^  ct  yx        "^oxy  y  ^^x        "  ox^  y 

CT^t.   -    Y   4.V      +  (Tuv      -   Y     uv      +   vv      =   -    ID      -  (T 

^    t        '^  ot   y  X        •'  ox     y  y  "7 

(2.11.)  (Tv^  -  y^^v^  -  ^y  =   0 

p  =   0     at     y  =  y^(x,t) 

J    .    +(Sy        +   uy        +CJ'uY.      =  V      at      y  =   Y=(x,t) 
''ot  -^ox  '-ox  ''IX  "^         .)  ^\     t     I 


The  para-aeter  (5*  in  equations  (2.'i)  will  be  small  if  the 
vertical  dimension  b  is  small  compared  to  the  horizontal 
dimension  a,  which  is  the  case  for  thin  jets.   In  such  cases 
we  ma^?"  seek  a  solution  of  equations  (2. Li.)  in  the  form  of  a 
power  series  in  (jr  ,      Tmis  we  let 

CO    ,   X  00    ,   > 

,(n),   „  ^\rrT^  ,,   -  K ,,(nj,^  „  ^\^n 


n=0  n^Q 

(2.5) 

CO    ,   V  00    ,   V 

P  =EIp^''^^>-»t)a'^  ,    y,  =ZZy ''^^,y,t)(cr"    . 

n=0  -^    n=0 


To  determine  the  coefficients  in  these  expansions  we  insert 
equations  (2,5)  into  equations  (2.3)  and  (2.[|.)  and  equate  the 
coefficients  of  each  power  of  (j   to  zero.   In  this  way  we 
obtain  from  the  coefficients  of  T  : 
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6. 
(2.6) 

a)  y^°'(x,t)  -   -y^°'(x,t) 

b)  -V  u^°)  +  v^°)  =  0 

■^ox  y      y 

N       (o)       (o)  (o)    (o)  (o)       (o) 

'     ot  y    '^  ox     y         y     "^  cx^  y 

.>       (o)      (o)  (o)  ^  (o)  (o)      (o) 

d)  -Yot^y  -  Yox^  ^y  ^^   'y    -Py 

e)  y  v^°)  -f  u(°)  -  0 

^ox  y      y 

f)  P^"^  =0     at  y  =  y:°^x,t) 

g)  Yot  ^  ^^°^yox  =  ^^°^  ^^  y  =  yi°^x,t). 

The  procedure  thus  far  is  nearly  the  same  .as  that  of  [2] 
and  [3j"   1'he  essential  differoncc  occiirs  in  the  choice  of 
dimonsionless  variables . 

3«   Zero- Order  Solution  for  Unsteady  Jets 

Prom  equations  (2,6b)  and  (2,6e)  it  follows  that 

(1.1)  u^^^  =  u^'^^x,t)   ,    v(°)  =  v^°^x,t)    . 

Now  using  equations  (3»1)  in  (2.6c),  we  have  p;^  '  =  0,  and 
then  from  (2.6f) 

(3.2)  p^*^)  =  0    . 

With  these  results,  equations  ( 2 .6b) -( 2,6f )  are  identically 
satisfied,  and  only  equations  (2.6a)  and  (2.6g)  remain. 
These  equations  are  Insufficient  to  complete  the  determination 
of  the  zero  order  terms,  and  therefore  we  proceed  to  equate 
to  zero  the  coefficients  of  <X      in  equations  (2,3)  and  (2.14.). 
Wo  obtain: 


.1....  -. 


'"\: 


7^ 
(3.3)  rx)  y[l\x,t)  =  -y^^^x,t) 


b)       u(^"^  -  V  u^l^  +  v^l)  =  0 
'       X     ■■'  ox   y      y 


0,  4°'.u(°)4°)  =  .„,p<.i' 


d)  vl°'  +  u'°'v'.°'  =  -pi^'-l 

L".  X  J 


e)      v(-^  -  y  v^l^  -  u(l^  =  0 
'       X     '^ox  y      y 


f)  p^"-  -  0     at  y  =  y^^^ 


.(1)  -  r,    „.  ..  -  ,-(o: 


,    (o)     (o)  (o)  ^   (1)       (1) 
g)  Jl^'    +  n'    'y|^'  +  u'  'y^^  =  v^  '   at   y  =  y 


(o) 

i 


( ■•  ) 

From  cauntions  (3.3c)  and  (3.1)  wo  find   that  p;,"^   is 

( -\\  J 

indciDcndent  of  y.   Since  p'  '  is  zero  at  two  different  values 

of  y,  by  equation  (3»3f)»  we  have 

(3.^)  P^"^'  =  0    . 

Then  equation  ['^*'}>o)    becomes 

(3.5)  u[°)  ^  u(°)u^^°)  =  0    . 

This  last  equation  implies  that  the  zero-order 
horizontal  velocity  component  is  a  constant  for  each  "particle" 
of  the  fluid.   This  suggests  the  following  method  for  obtaining 

u^'^^  Lot 

(3.6)  u^°)  =  A(r) 

where  A  is  an  arbitrary  function  of  a  parameter  t!  .  This 
parameter  ^  is  defined  implicitly  in  terms  of  x  and  t  by 
the  equation 

(3.7)  X  =  (t-tr)  A(-2:)  +  x'"*(t^) 


?/).. 
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Here  x'"('C)  is  also  an  arbitrary  function.   We  may  interpret 

V    as  the  time  at  which  the  particle  at  x  and  time  t  was  at  x""''. 

For  example,  x'^v  C)  might  be  the  x-coordinate  of  the  nozzle 

at  time  V  ,    and  the  emitted  particles  would  then  have  the 

vjjlocity  A('t;).   It  is  easily  verified  by  differentiations 

that  equation  (3 •6)  gives  a  solution  of  equation  (3.5) > 

provided  t?  Is  defined  by  equation  (3.7). 

In  the  same  way,  a  solution  of  equation  (3.3d)  for 
( o) ,   .  %  . 

V  ^   '  ( X , t )   IS 

(3.8)  v^°^  =  -{t-Z)  +  w(tr) 

where  "^'liz:)     is  an  arbitrary  function.   Now  we  may  solve 
equation  (2.6g)  for  y  (x,t)  by  assuming  that  y  (x,t)  =  B(-2^,t) 
i.e.  y^  is  a  function  of  V   and  t.   Then  equation  (2,6g) 
becomes 

(3.9)  B^  =  v(°^' 

Using  equation   (3.8)    in  cqtxation   (3.9)    and   integrating  vie 

have    for  B   or  v 
"  o 

(3.10)  yQ(x,t)  =  y-=^(t^)  +  (t-r)¥(f )  -  ^^"1^ 

Here  y"""(  "K)  is  an  arbitrary-  function  vfhich  can  be  interpreted 
as  the  y-coordinate  of  the  nozzle  at  time  ^  . 

Before  considering  y\^°^(x,t),  v;o  first  integrate 

equation  (3.3b)  from  y|°'  to  y^°'  and  obtain: 


(3.11)       u('^y^°^  -  ■A'^h    =u(l)y^^  -  v^l) 


1„(o) 
(o) 

V 

^1 


Now,  subtracting  equation  (3.3g)  o-t  y  =  y^^  from  the  same 

(  o  i  ■'• 

equation  at  y^  ',  and  employing  equations  (3.3a)  and  (3.11), 

wc  have 

(3.12)  y(^)  .  (u(°-'y|°))^  =  0    . 
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This  equation  for  y^    can  also  be  solved  by  assioining 

y^  '  (z,  t)  =  P('£r,  t)  .   We  have,  after  making  this  substitution; 

(3.13)  P^  +  p4°^  =  0    . 

This  is  a.  first  order  ordinary  differential  equation  and  can 
be  solved  Immediately  yielding  for  P  or  y-j   : 


(3.11;)      y[°\x,t)  =  -y(°^x,t)  =  -—^ 


x'^+{t-  'C)A  -  A 


Here  C('C)  is  an  arbitrary  function,  which  may  be  related  to 
the  position  of  the  wall  of  the  nozzle  (or  to  the  nozzle  width) 
at  time  '^  . 

Vie  have  now  determined  all  the  zero-order  terms  (and  one 
first  order  term)  In  the  solution.   Collecting  our  results, 
we  have  in  zero-order: 

u  =  A(2r)  +  ... 

V  =  -(t-r)  +  ii(r)  +  ... 

(3.15) 

p  =  0  +  CT-Q  +    ... 

y^  =  7"''"(r)  +  (t-2^)  w(r)  -  i^::Jl- 

Jl  =  -72  =  C(r)  [x'^+  (t-t:')A^-A]~^ 

Xifhere    "C    is   defined  by 

X   =    (t-^)  A(tr)  +   x""'(tr) 

The  fijinctions  x''%  y''%  C,  A,  VJ  are  arbitrary.   If  the 
additional  conditions  on  the  flow  are  that  it  emerge  from  a 
nozzle,  then  x"  and  y''"'  must  be  taken  as  the  coordinates  of  the 
center  of  the  nozzle,  and  C,  A,  W  are  chosen  to  fit  the  width 
and  velocity  of  the  jet  at  the  nozzle,  all  at  time  -£r  .   We 
see  that,  to  this  order,  only  the  average  velocity  over  the 
cross- section  can  be  imposed  as  the  boundary  condition.   This 
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Is  due  to  the  asymptotic  character-  cf  the  approximation.   In 
the  higher  order  solution  iriore    details  of  the  velocity  at  a 
cross-section  can  bo  imposed.   We  must  also  point  out  that 
the  solution  fails  at  points  where  y  (x,t)  intersects  itself. 

The  solution  (3»l5)  resembles  the  flow  of  independent 
particles,  each  of  which  maintains  a  constant  horizontal 
velocity  and  describes  a  parabolic  trajectory.   The  pressure 
is  zero  (even  to  first  order)  but  the  width  of  the  jot  varies, 
showing  a  "bunching"  effect.   For  example,  the  jet  thickness 
(for  fixed  t^   )  increases,  remains  constant,  or  decreases 
according  as  A^  is  positive,  zero,  or  negative.   If  gravity 
is  neglected,  the  only  change  is  the  omission  of  the  last  term 
(one)  in  equation  (3.3d)  and  consequently,  the  solution 
equations  (3.15)  still  holds  provided  that  we  omit  the  first 
term  in  the  expression  for  v  and  the  last  term  in  the 
expression  for  j    . 

If  the  flow  is  steady  (Independent  of  t)  x",  y",  C,  A, 
and  W  are  constant.   Then  the  solution  equations  (3.15)  holds 
provided  that  vje  replace  t-  'V  by  xA~   from  equation  (3»7). 
The  resiAlt  is 

u   ~   A    +    ... 

V  =  W  -  xA    +  . . . 


(3.16)  p  =  0  +  (fc  +  . .  , 

•/q  =  y"  +  WA  X   -  ^  A    -X 
y^  =  -y2  =  CA"   +  . . . 


In  this  case  the  vertical  width  of  the  jet  remains  constant 
to  this  order.   This  solution  does  not  apply  for  vertical  jets, 
in  which  case  A  =  0,   A  somewhat  different  method  must  be  used 
to  treat  this  case,  and  it  is  described  in  Appendix  I. 

The  zero  order  solutions  obtained  in  this  section  are 
exactly  the  same  as  the  results  of  the  hydraulic  theory.   Thus 
the  hydraulic  theory  yields  a  first  approximation  to  the 
solution  of  the  hydrodynamic  problem,  valid  for  thin  jets. 
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4-'   First  Order  Solutionf  or  Steady  Jets 

hie  shall  novj  proceed  to  the  detorminatlon  of  the-  first 
order  terms,  but  because  of  the  complexity  of  the  equations, 
v/e  shall  consider  only  the  steady  state  solution. 

From  equations  (3.3t)»  <^   and  g)  one  easily  finds  that 

1/        \  1/        0\  <^/n   .   /   2   ^2^-l,  Ox 

u  (x,y)  =  u  (x,y^)  +  v^(l+  {p,^^)  )   (y  -y-^; 

(I^.l) 

v-'-(x,y)  =  y^,^u  (x,y) 

Equations  (l+.l)  and  the  preceding  expressions  for  p   and  the 

zoro  order  solutions  satisfy  equations  (3»3b-g).   To  complete 

11  1 

the  determination  of  u  and  v  and  to  find  y  ,  it  is  necessary 

2 

to  equate  to  zero  coefficients  of  C   in  equations  (2.3)  and 

(2J4.).   From  the  second  and  third  equations  resulting  from 

equations  {2.1\.)    we  eliminate  pA.   and,  by  substitutions  from 

equation   (i|-.l)  we  obtain  an  ordinary  differential  equation 

in,  u  .   This  equation  is  easily  solved  and  thus  (from  equation 

(i|..l)  )  both  u"  and  v  are  found.   Now  employing  these  results 

in  the  eqiiations  resulting  from  the  first  and  sixth  of 

equations  (2.L|.),  we  obtain  y.  .   This  completes  the  determination 

of   the   first   order    solvition. 

The  resulting  solution,  including  both  zero  and  first 

order  term,  is: 

u  =  A+O'i  A^-A^y )  ( A^  +  a\'^  -  2AWx  +  x^ )  ""^  +  .  .  . 
V   =  W-  A"-'-x  +  cr(A^D-Ay)(A^"^+ A^W^-   2AWx  +  x^ )  "  "'■  ( AW   -x)  +   ... 
()|.2)      p   =   0  +  (T'O  +... 


■^2 


=   -CA'-'-    -  (rCA^D(  A^  +  aSj^  -  2AWx  +  x^  )"^  + 


Here  D  is  a  constant,  which  like  the  other  constants,  may  be 
determined  from  the  prescribed  boundary  values  of  the  solution 
at  some  cross- section  of  the  stream. 
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5 .   Unsteady  Flow  on  a  Rlgi.d  Surface  (o.g.  Hill  or  Spillway ) 

In  the  preceding  problr;m  the  jet  had  two  free  stream- 
lines.  But  in  the  flow  over  a  rigid  surface,  such  as  a  hill, 
spillway,  river  bed,  etc.  there  is  only  one  free  streamline 
(except  when  the  flow  detaches  itself  from  the  surface).   Such 
flows  with  one  free  streamline  can  be  treated  by  the  ^receding 
method,  with  a  slight  modification.   In  fact,  it  was  for  the 
treatment  of  a  problem  with  one  free  streamline  that  the 
method  was  originally  employed  -  namoljr  th<j  problem  of  waves 
in  shallow  water.   Hov:cver,  in  that  problem  the  "bottom"  was 
ass\imcd  to  be  practically  horizontal,  and  its  slope  small. 
In  the  present  nroblom  we  Xi?ish  to  consider  arbitrary  "bottom" 
surfaces  which  may  have  largo  slopes.   For  such  cases  the 
original  discussion  of  the  shallow  water  theory  does  not  apply. 

We  formulate  the  problem,  just  as  in  section  2.   However, 
the  last  two  of  oquations  (2.1)  apply  only  at  the  single  free 
surface  y  =  y  (x,t).   At  the  "bottom"  surface  y  =  ^n^^^  which 
is  assumed  to  be  known,  we  have  the  condition 

(5.1)  y_u=v  at   y=  y^(x) 

ox 

The  dimensionlcss  variables  in  equations  (2.2)  and  the 
equations  (2.Li.)  arc  the  same.   Hovjever,  jA^)    =  ay^(x) 
represents  the  bottom  surface  instead  of  the  center  line  as 
in  the  previous  case.   The  derivative  y  .^  is  also  identically 

O  u 

zero.   Equation  (5«1)  becomes 

(1^.2)  uy    =  V   at   y  =  0    . 

The  scries  representation  of  the  solution,  equation  (2.5)  and 
the  equations  (2.6)  are  the  same  as  before.   The  method  of 
solution  is  also  similar  to  that  of  the  previous  case. 

The  zero- order  solution  (with  p  given  to  first  order) 
is  found  to  be; 
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ti 

)  -2yol 

_  1 

+ 

2 

11/2 


V 


k{^)  -2y 


-2v^V2 


o 


L  1  +  y 


T 


ox 


(5.3)     P  =  0+cr(y5^°--y)(l  +  y^^+y^^^[k(-2;)  -2y^] )  (1+ Yo.^)"^  +•.. 


7, 


n^)z 


X 


[  1*4  1 

+   ...   *=  i :  V   t.  j  I  — ; — - — - — 


1/2 


■1  +  -^  (1  + 


^ox^^^^^^^'^^-^^^o)"^'^^'^^ 


+    X. 


2    ,    -;>^>      -il/2 
1  +  Y      ( X    )  ' 

"^  ox  ^       ' 


Lk(tr)  -2yQ(x-'')J 


>  + 


where     ^C   is   ciefined   as   a    function   of  x   and   t   by 


(5.1^) 


t  -  r  = 


X 


x-"(r) 


k(  f)  -2y^{a) 


1-1/2 


L  1  +y      (a) 

-^  ox  ^    '     -^ 


da 


Here  k('C),  x'^CK)  and  ¥,{t:)     are  arbitrary  functions. 
The  flow  resembles  a  stream  of  individual  particles  sliding 
along  the  rigid  surface.   Then  V:{li)    and  x"{'2:')  represent  the 
kinetic  energy  and  x-coordinate  at  time  "tr  of  the  particle 
which  is  at  x  at  time  t.   The  depth  y-,  of  the  stream  exhibits 
a  "bunching"  effect  as  in  the  previous  problem.   Even  to  this 
order  mass  and  energy  are  conserved. 

The  pressure,  given  by  the  third  of  equations  (5.3) j  is 
approximately  hydrostatic  where  the  bottom  slope  y   is  small. 
Hov;ever,  if  the  slope  is  large,  the  variation  of  pressure  v;ith 
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depth  is  much,  less  than  vhat  would  be  given  bj   hydrostatics. 
If  the  pressure  at  somo  point  on  the  bottom  (y  =  0)  becomes 
zero,  (and  would  be  negative  for  neighboring  points)  then  the 
streari  detaches  itself  from  the  bottom  at  this  point,  and 
becomes  a  stream  x^7ith  two  free  streamlines,  to  which  the 
previous  solution  applies.   This  will  occur  if 

(5.5)  1+V^  +y    (k(^)-2v  )  =  0 
V  -■  ^  /  -^^  ox  "  oxx  ^       "  o 

Thus  the  stream  detaches  vjhen  the  bottom  curvature 
becomes  too  great.   In  fa.ct  the  curvature  at  the  point  of 
detachraent  is  (in  the  original  variables) 

(5.6)  y  __(l  +  y2_)-3/2  =  _g--2^^^ -2^^3/2    ^ 

oxx     ox  ox 

This  is  the  curvature  of  tho  trajectory;'  of  a  freely  falling 

particle  with  horizontal  velocity  u  and  trajectory  slope  y 

ox 
Thus  we  see  that  if  the  curvature  of  the  bottom  exceeds  that 

of  the  path  of  a  freely  falling  particle  with  the  same 

horizontal  velocity  and  slope  as  the  stream,  detachment  occurs, 

(If  the  pressure  is  zero  but  not  negative  on  the  bottom,  it 

is  net  necessary  to  speak  of  detaclimcnt  since  the  "detached'' 

stream  would  also  have  the  bottom  as  its  streamline.) 

If  tho  flow  is  steady  ( independent  of  t)  the  above 

solution  (5.3)  is  valid  with  k,  x'"'  and  IT  constants.   Equation 

{^,l\.)    is  then  unnecessary,  and  the  remaining  discussion 

applies . 

6.   The  Effect  of  Surface  Tension  on  a  Jet 

In  order  to  take  account  of  surface  tension  T,  we  need 
merely  modify  the  fifth  of  (2.1),  which  now  becomes 

(2.1)'"'      ?=(-l)^TK^    at  y  =  y^(x,t)  ,   (i  =  l,2)   . 

The  sign  is  determined  by  requiring  that  i  =  1  corresponds  to 
the  upper  surface,  and  K.    is  the  curvature  of  surface  i.   In 
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the  introduction  of  dimensionless  variables  we  include  T,  K. 
defined  by 

T  =  pgabT   ,    K.  =  a  "^K. 
Then  the  fifth  of  (2.1^.)  becomes 

(2.I4.)""'         p  =  (-l)^(rTK^   at   y  =  y^(x,t) 

The  subsequent  expressions  given  by  {2,$)    and  the  zero 
order  equations  (2.6)  and  their  consequences  (3»1))  (3«2) 
remain  unchanged.   The  first  change  occurs  in  {3"})    f)>  which 
becomes 

(3.3)  f)'^      P^^^  =  (-1)^  Tk(°)    at   y  =  y(^^    . 

As  before  (3.3)  c)  and  (3.1)  Imply  that  v[r        ^^ 
Independent  of  y.   Corribining  this  fact  with  (3.3)  f)   above, 
we  find 

{^..lir  P^^^  =  -TK^°)  -^^-v    . 

^1  ' 

This  replaces  the  former  (3.)^).   Using  (3,[|.)",  (3.3)  c)  and 
(3 •3)  <^)    become 

(3.5)"  u4°^  +  u^^°^u^°^  =   -  TK^°^  -if 

t  X         (  o) 

■^1 

These  two  equations,  along  with  (2,6)  g)  which  is 

(2.6)  g)  y^^  +  ^^-^"^ox  =  ^^°^ 

provide  three  equations  for  the  determination  of  u^   ,  v^°  , 
y^  and  j\° ' •      The  fourth  equation  is  (3.12),  as  before 
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(3.12)  At     "-  ^''^°^A"^^^  =  ^       • 

I^ewriting  (3.5)  ,    {3»6)  ,  (2.6)  g)  and  (3.12),  omitting 
all  indices  and  letting  h  =  yi  'we  have 

Tv 
v^  +  uv   -    "^^ 


t      X    ,  ,,   2.3/2 
'  ^x' 


yt  ^  ^y^  ^  ^ 

h^  +  (^*)x  =  0   . 

This  is  a  set  of  four  partial  differential  equations  for  the 
four  functions  u,  v,  y  and  h,  which  depend  upon  x  and  t. 
Exactly  the  same  set  of  equations  is  obtained  in  the  hydraulic 
theory,  where  the  jot  is  treated  as  a  stream  of  individual 
particles  with  mass  proportional  to  the  jet  v;idth,  provided 
the  force  of  surface  tension  is  included  as  well  as  gravity. 
In  that  case  however,  the  left  hand  sides  may  be  thought  of 
as  total  time  derivatives. 

'7e  will,  now  solve  these  equations  when  the  jet  is 
steady,  i.e.  independent  of  time.   In  this  case  the  last 
equation  above  can  be  integrated,  which  introduces  the 
constant  flux  c.   Then  the  remaining  equations  can  be 
simplified,  yielding 

uh  =  c 

V  =  uy^ 


1     / ,. 


4;' 


(■■-^^..    '> 


I.'  V'^  ^ 
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The  third  of  these  eqiiations  can  be  Integrated  from  0 
to  X,  and  yields 

c(u-u{c))  =  T[(l+y^)"^/^-  (l+y^(o))"^/^] 

Using  this  to  eliminate  u  from  the  last  preceding  equation, 
we  obtain  a  second  order  equation  for  y.   Introducing 

a  -   cu(o)  -  T(l+yf(o))"  '       this  becomes 


[T(l+y^)-^/2  +  a]  =  -c^    . 


^■■^xx '  ~  '  ~  "  X 

Integrating  from  0  to  x  \^&   have 

-1-12  2 

aT[sinh  y^  -  slnh  7^^(0)1  +  a  [y^  -  y^(o)]  =  -c  x 

The  last  equation  is  a  first  order  equation  for  y(x) 
involving  several  para;i':eters .   It  seems  that  this  equation 
must  be  integrated  nunaerically,  and  it  is  therefore  important 
to  reduce  the  number  of  parameters.   To  this  end  we  may  set 
y  (o)  =  0,  since  the  resulting  curve  will  apply  to  a.rc<r   Initial 

slope  if  read  from  that  slope  on.   Next,  we  introduce  the  new 

2  2 

c  c 

variables  — ^  x  =  £  and  — p-y  =  ">^  .   In  terms  of  these  variables 

a  a 

the  equation  becomes 

\   +  \   slnh-1  ^^  =  -^    . 

This  equation  involves  only  one  param.eter,  —  =  [— ~-^-l]~  , 

and  is  therefore  convenient  for  nuraerical  integration.   V/hen 

p 
T  f 

—  =  0,  the  aoove  equation  yields  the  expected  parabola  V3  =  - '^ 

In  figure  1  graphs  of  the  solution  of  the  above  eauatlon  are 

T 
shown  for  various  values  of  —  corresponding  to  various  values 
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cuf  o )  '^ 

— 1^:; — -  .   In  the  E,,  yi    ^-ariables  the  curves  vjith  —  >  0  lie 

T 
above  the  parabola  obtained  for  —  =  0.   However  in  the  y.,j 

variables  these  curves  lie  below  the  parabola.   Since  this 

is  the  physical  plane,  it  is  important  to  recognize  that  in 

the  physical  plane  the  jets  with  surface  tension  lie  below 

T'  GUI  o) 

the  parabola  provided  —  >  0,  or  ~~ — -  >   1. 

T 
It  is  interesting  to  notice  that  —  becomes  negative 

CU(o)   ^  ,     ,        ,      4.'     1  •    J.  •  1    -n   2,   > 

w^ien  — ^ — -  <   1,  I.e.  when  tne  kinetic  energy  ^»  2hu  (o) 
becomes  less  than  the  surface  energy  T.   In  this  case  the 
above  equation  yields  the  result  that  the  jet,  although 
starting  horizontally,  vjill  at  first  rise  until  it  reaches 
a  point  at  which  the  curvature  becomes  infinite  and  the  theory 
consequentl;/  fails  (see  Figure  1).   This  surprising  phenomenon 
can  occur  only  in  extremely  thin  and  slow  jets  since  the 
surface  tension  is  very  small  for  all  ordinary  liquids.   In 
such  slovj  jets,  near  the  orifice  or  nozzle,  the  hydraulic 
theory  may  not  apply,  and  therefore;  the  effect  night  not 
occur.   However,  in  a  fluid  with  a  large  surface  tension, 
i-:''  which  tais  effect  could  be  expected  for  a  reasonable  jet 
velocity  and  thickness,  Lhc  hydraulic  theory  should  apply 
and  the  effect  should  occurs.  Kaen  ^-^j—-  =   1,    the  steady 
state  equations  have  no  solution,  and  therefore  the  hydraulic 
theory  cannot  apply  in  this  special  case. 
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Appendix  I 

yertic:-il.  Steady  Jets 

We  now  consider  a  jet  which  flows  vertically,  and  for 
simplicity,  is  steady.   The  method  of  sections  2-l,\.   is  not 
applicable  in  this  case,  and  therefore  we  will  describe  an 
alternative  method.   The  problem  is  formulated  just  as  in 
section  2,  but  the  free  streamlines  arc  designated  by 
X  -   x.(y),  i  =  1,2.   Equations  (2.1)  all  remain  valid  except 
the  last  which  becom.es 

(A.l)         Vx    =  u   at   X  =  x,(y)    (i  =  1,2)   . 

A   suitable    transformation  to   diraensionless   variables   is 
given  by 

X  =  ax  u  =  Ssh  u  X.    =  ax, 

(A. 2) 

y   =   by  V   =  Tgb  v  p   =   gbpp 

I-Jo   now  proceed  as    in   the   previous    case,    but   v;ith    CT  =  •^  . 

'  b 

The  result  of  determining  the  zero  and  first  order  terms,  for 
a  symmetrical  let  (x-,  -   -Xp )  is: 

u  =  0-  Gr(k^-2y)"^/^x  +  ... 

V  =  -(k^-2y)^/2_^^^^^^j^_^_2y)-l/2  ^  ^^^ 

(A. 3) 

p  =  0+G'.0+CJ^[k^(k^-2y)~2- x^(k-^-2y)~^]  +  ... 

x^  =  -^2  "  -k2(k-L-2y)"^--^2+<y[-k|^(k^-2y)~^-/2+k2k^(k^_2y)-3/2j  +  _, 
Here  k^ ,  kp,  k^,  k.  are  constants. 
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By  Introducing  the  original  variables,  we  obtain 
u  =  -xi(k^-2y)~^/^+  ... 

(A.[|.) 

^1  "  -"^2   =  -{k2+k^)(k;^-2y)"^/2  +  k2k2(k^-2y)'^/2^  _^ 

2     -  2      -^     -  -1 

P  =  pg  ^^(k^-gy)   -  pgx^(k^-2y) 

In  equations  (A.I4-),  k^-k,  are  constants  to  be  determined 
frora  the  boundary  conditions  at  some  cross- section  of  the 
.jet.   The  last  form  of  the  solution  was  given  to  show  that 
the  particular  choice  of  a  and  b,  or  O*  ,  does  not  affect  the 
solution.   This  is  also  the  case  with  our  previous  resv'lts. 

The  jet  is  symmetric  about  the  vertical  center  line 
and  its  vertical  velocity  varies  approximately  like  that  of 
a  freely  falling  body.   Its  v/idth  varies  inversely  as  the 
vertical  velocity,  as  we  would  expect  from  conservation  of 
mass.   A  small  horizontal  velocity  toward  or  away  from  the 
center  of  the  jet  acco"ants  for  its  thinning  (when  falling) 
or  broadening  (when  rising) .  (The  positive  square  root  applies 
to  falling  jets;  the  negative  square  root  must  be  used  for 
rising  jets.) 
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